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An effective quasiparticle description of the thermody- 
namics of deconfmed matter, compatible with both finite- 
temperature lattice data and the perturbative limit, is gener- 
alized to finite chemical potential. Within this approach, the 
available 4-flavor lattice equation of state is extended to fi- 
nite baryon density, and implications for cold, charge-neutral 
deconfmed matter in /3-equilibrium in compact stars are con- 
sidered. 

Key Words: quark-gluon plasma, deconfinement, equation of 
state, quasiparticle 
PACS number(s): 12.38-Mh, 31.15.Lc 

The equation of state (EoS) represents an important 
interrelation of state variables describing matter in local 
thermodynamical equilibrium. All microscopic charac- 
teristics are integrated out and only the macroscopic re- 
sponse to changes of the state variables as e. g., the tem- 
perature T and the chemical potential \x are retained. 
Via the Gibbs equation, e = Ts + fm — p, the energy 
density e is related to the entropy and particle densities, 
s = dp/dT and n = dp/d/i, respectively, and the pressure 
p. The thermodynamical potential p(T, fj) thus provides 
all information needed to evaluate, e. g., sequences of stel- 
lar equilibrium configurations by means of the Tolman- 
Oppenheimer-Volkov equation, or the evolution of the 
universe via Friedman's equations, or the dynamics of 
heavy-ion collisions within the framework of relativistic 
Euler equations. In the examples mentioned not only ex- 
cited hadron matter is of relevance, rather, at sufficiently 
high density or temperature, a plasma state of deconfmed 
quarks and gluons is the central issue. 

Quantum Chromodynamics (QCD) is nowadays the 
generally accepted fundamental theory of interacting 
quarks and gluons. The challenge, therefore, consists in 
the derivation of the EoS of deconfmed matter directly 
from QCD. In first attempts JLj the EoS of cold quark 
matter was derived perturbatively up to order 0(g 4 ) in 
the coupling g. Finite temperatures have been consid- 
ered [|2|, where the perturbative expansion is extended 
up to the order 0(g 5 ). However, in the physically rele- 
vant region the coupling is large, so perturbative methods 
seem basically to fail and, consequently, nonperturbative 
evaluations are needed. Present lattice QCD calculations 
can accomplish this task, and indeed the EoS of the pure 
gluon plasma H and of systems containing two or four 
light dynamical quark flavors HJq] are known at finite 
temperature, and simulations with physical quark masses 
may be available in the near future. 



As a matter of fact, current QCD lattice calculations 
are restricted to the chemical potential fi — while a 
more detailed understanding of, e.g., the structure of 
quark cores in massive neutron stars, the baryon con- 
trast prior to cosmic confinement, or the evolution of 
the baryon charge in the midrapidity region of central 
heavy-ion collisions requires the EoS at vanishing net 
quark density. Here we are going to develop an EoS, 
based on a phenomenological quasiparticle model, which 
extrapolates available lattice results at \x = to finite 
values of /i (even at T = 0) and at the same time in- 
terpolates smoothly to the asymptotic regime of QCD 
at T — > oo and \x — * oo, with the strangeness degree 
of freedom properly included. With our model we focus 
on deconfmed matter and do not touch upon problems 
of confined matter or the deconfinement transition itself; 
the latter might be covered via usual ad hoc procedures. 

We consider an SU(N C ) plasma of gluons, with N c — 3 
for QCD, and Nf quark flavors in thermodynamic equi- 
librium. Within the approach outlined below, the inter- 
acting plasma is described in terms of a quasiparticle sys- 
tem, a picture arising asymptotically from the in-medium 
properties of the constituents of the plasma. For ther- 
mal momenta k > 0(T) the relevant modes, transversal 
gluons and quark particle excitations, propagate predom- 
inantly on-shell with dispersion relations u>f(k) rs mf+k 2 
and 
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while the longitudinal gluonic and helicity flipped quark 
states are essentially unpopulated ||. Neglecting sub- 
leading effects, n* are the leading order on-shell selfen- 
ergies of parton species i. Depending on the coupling 
G 2 , the temperature and chemical potential as well as 
the rest mass mot (the latter vanishing for gluons), the 
n* are given by the asymptotic values of the hard ther- 
mal/density loop selfenergies B , 



n* = 2 UJ qQ (m + UJ q o) , UJq Q 
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Generalizing the approach of ref. J/J to a finite chemical 
potential /i controlling a conserved particle number, the 
pressure of the system can be decomposed into the con- 
tributions pj of the quasiparticles, and their mean field 
interaction B, 



p(T, u- ml,) 



v^T^^-m^-BCn*) 



(3) 



where formally pj is the pressure of an ideal gas of bosons 
or fermions with state dependent effective masses (|[||). 
By the stationarity of the thermodynamic potential Q = 
—pV under functional variation of the selfenergies ||, 
which in the present approach simplifies to dp/dH* — 0, 
B is related to the quasiparticle masses, 



dB dp 3 (T, Hj ; m 
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which implies that the entropy and particle densities are 
given by the sum of the quasiparticle contributions 



dp j (T, f i J ;m 2 ) 



dT 



dp J (T,ii j ;m 2 j ) 



dfij 
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The quasiparticle approach represents an effective resum- 
mation of the leading-order thermal contributions ||. 
Hence, it is expected to be an appropriate framework as 
long as the spectral properties of the relevant plasma ex- 
citations do not differ qualitatively from their asymptotic 
form. In the hot </> 4 theory, for which an equivalent ther- 
modynamic quasiparticle description was derived by re- 
summing all tadpole diagrams En] , this assumption was 
shown to be fulfilled also at larger values of the coupling 
by resumming the propagator beyond 1-loop order El| . 
Corroborating the assumption for QCD thermodynam- 
ics at li — 0, the quasiparticle approach was found to be 
in a remarkably good agreement with lattice data even 
close to the confinement region |12-|l4] , with the effective 
coupling G 2 in eq. (En parameterized by 



G 2 (T, M = 0) 
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interpolating smoothly to the asymptotic QCD limit. 

For finite values of ii, we observe that due to the sta- 
tionary property (||) the dependence of the function B 
on the state variables is determined by 
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The pressure (||) and thus the function B being a poten- 
tial of the state variables T and li, the functions Bt and 
B^ have to respect the integrability condition 
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In the selfenergies, eq. (||), the effective coupling itself is 
a function of the state variables, so eq. (0) represents a 



first order partial differential equation for G 2 (T, u) . With 
G 2 {T,li) given, e.g., at n = by lattice data, this flow 
equation determines the effective coupling and hence by 
eqs. (|l|-||) and (|7]) the EoS of the plasma at finite tem- 
perature and chemical potential. 

As a first example let us now apply the quasiparticle 
approach to the EoS of the QCD plasma with Nj = 4 
light flavors, which is numerically known for a — in a 
restricted interval of T. From a fit (cf. fig. 1 in ref. E3[) of 
the lattice data [[| we obtain the parameters A = 6.7 and 
T s /T c = —0.81 of the effective coupling ([j|), and the gluon 
and quark quasiparticle degrees of freedom d g = 20.6 



and d„ 
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d g . The large value of the parameter 



"1 ~ 2(JV C 2 -1) "»• 

T s /T c points to the distinct nonperturbative behavior at 
T — > T c . We note that the ratio of the fitted value of 
d g to the expected value of 16 is of the same order of 
magnitude as finite lattice size effects in the interaction 
free limit. 

With this at hand we extrapolate the lattice data to 
finite values of ii q = /^, corresponding to a finite net 
quark density. Eq. flq) is a quasilinear partial differential 
equation of the form 
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with the coefficients ar^ and b depending on T, li and 
G 2 , which can be solved by the method of characteristics. 
It is instructive to consider first the asymptotic limit, 
m 2 /T 2 ~ G 2 — > g 2 — > 0, where (||) reduces to 
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= 0, (10) 



with c = (AN C + 5N f )/(9N f ). This equation yields 
g 2 = const along the characteristics given by cT A + 
2T 2 (li/it) 2 + (/i/ir) 4 = const. As a result of this elliptic 
flow, the renormalization scale T c /X of the effective cou- 
pling g 2 (T, 0) determines the scale T c 7rc 1 / 4 /A of g 2 (0, li), 
and the pressure (0-0) coincides with the perturbative 
QCD expression up to the order 0{g 2 ). 

The flow (g) of the effective coupling, with G 2 (T, 0) ob- 
tained from the lattice data referred to above, is shown 
in fig. yj. The characteristics related to larger tempera- 
tures resemble their asymptotic form, and the coupling 
displayed in fig. || decreases asymptotically as expected. 
For smaller values of T and a, the pronounced increase 
of G 2 is considered to be a sign of the vicinity of the 
phase transition. Indicated by intersecting characteris- 
tics, the solution of the flow equation is non-unique in a 
certain low-temperature region. In this region, however, 
the pressure turns out to be negative, so the ambiguity is 
of no physical relevance. To obtain the pressure, which 
is available in tabular form frlql , the function B(T,u) 
can be computed along the characteristics using eq. (H). 
Hereby, the required function B(T,u = 0) is determined 
by the first equation of (pj) up to an integration constant 



Bq = B(T C , 0) which is fixed by equating the quasiparti- 
cle pressure and the lattice data at T c . For chemical po- 
tentials [i ~ 2.5 T c and small temperatures, in a region of 
the phase space corresponding to characteristics emanat- 
ing from the vicinity of the point (T c , ix = 0), the pressure 
becomes negative. This region of instability provides a 
lower boundary for the value of the confining chemical po- 
tential fi c > 3.4 T c at T = 0. As evident from fig. [j], this 
boundary also excludes the region of non-unique flow. 
Hence, the quasiparticle model is intrinsically consistent 
and therefore considered to provide a realistic extension 
of the EoS obtained by lattice calculations to /i ^ even 
near the confinement transition. 

To conclude this example we remark that by the 
Fcynman-Hcllmann relation the quasiparticle pressure 
(0) leads to a chiral condensate (ijjip) ~ mo q , so for mass- 
less quark flavors the restoration of chiral symmetry in 
the deconfined phase is inherent in the approach. 

Thermodynamic lattice simulations of QCD with phys- 
ical quark current masses, too«,c2 ~ for the light fla- 
vors and m 0s ~ 150 MeV for strange quarks are still 
lacking, so the parameters of the quasiparticle descrip- 
tion cannot be fixed at finite temperature to extend the 
physically relevant EoS to \x ^ 0. However, a trial EoS 
p trial (T, \x = 0) can be constructed by reasonable choices 
of the model parameters until they may be specified by 
precise lattice data. Due to the stationarity property (0) 
of the thermodynamic potential, p trial (T, \i = 0) is ex- 
pected to possess only a weak dependence on the param- 
eters. The model parameters are restricted to match the 
hadronic pressure p had at the confinement temperature 
which we assume to be T c = 150 MeV. The uncertainty 
of p had (T c , 0) = 3.1 • 10 8 MeV 4 as obtained by a hadron 
resonance gas model can be absorbed into the integra- 
tion constant Bq, which we vary independently beside 
the scaling parameter A in eq. (||). Being related to the 
QCD scale parameter Aqcd, A is expected to be a slowly 
increasing function of the number of flavors; considering 
3 < A < 9 we cover the values A Ar/=0 = 4.2 as obtained 
from the pure gluon lattice data || and \ Nf=4 = 6.7 in 
the example above. 

With the aim to study implications for quark matter 
stars, we consider in the following a plasma of gluons, 
quarks and electrons in /3-equilibrium maintained by the 
reactions d, s <-► u + e + v e , which imply the relations 
fid = [x s = (X u + fx e = fx among the chemical potentials. 
The electron chemical potential /x e as a function of T and 
ix is determined by the requirement of electrical charge 
neutrality. With the electron pressure p e approximated 
by its free limit, the pressure p — p tIlal + p e and the re- 
sulting energy density are shown for several values of the 
parameters A and Bq at /x = in the left panels of fig. |3J. 
Already at temperatures slightly above T c , the scaled en- 
ergy density reaches a saturation-like behavior at some 
90% of the asymptotic limit. This feature, known qual- 
itatively from the lattice simulations [§-g, is to a large 
extent insensitive to the specific choice of the free param- 
eters, while Bq has a distinct impact on the latent heat. 



At vanishing temperature, the resulting EoS is displayed 
in the right panels of fig. |[ In this case, the asymptotic 
values are approached more slowly due to the less rapid 
decrease of the effective coupling with increasing values 
of fx, similar to the Nf — 4 plasma. 

With the resulting EoS e(p), sequences of hydrostatic 
equilibrium configurations of cold stars can be calculated 
by the Tolman-Oppenheimer-Volkov (TOV) equation (cf. 
eq. (2.212) of ref. |lq]). For energy densities up to sev- 
eral times the nuclear density and at temperatures less 
then some 10 MeV, the EoS of /3-stable quark matter, as 
estimated by our quasiparticle approach can be parame- 
terized by e = AB + ap. While the naive bag model EoS 
has a = 3, our quasiparticle EoS, with the considered 
choices of the model parameters, yields 3.1 < a < 4.5, 
indicating the nontrivial nature of the interaction. For 
parameters A > 5, values of B 1 ' 4 > 200 MeV are found, 
while only the extreme choice of A = 3 allows B 4 ! 4 as 
small as 180 MeV. However, the value of B, i. e., the en- 
ergy density at small pressure, has a strong impact on the 
star's mass and radius, obtained by integrating the TOV 
equation. Compact quark stars of mass M > M©, e.g., 
could only exist for values of B 1 / 4 < 180. ..200 MeV, de- 
pending marginally on a. We hence conclude that com- 
pact stars with typical masses ~ M© and radii ~ 10 km 
are unlikely to be composed purely of deconfined matter 
in /3-equilibrium, irrespective of the uncertainty of the 
model parameters of our approach. 

In summary we have generalized a thermodynamic 
quasiparticle description of deconfined matter to finite 
chemical potential fx not accessible by present lattice cal- 
culations. Of central importance to the model is the ef- 
fective coupling G 2 (T, fx) which can be obtained at ix = 
from available lattice data, proving at the same time the 
applicability of the effective description even close to the 
confinement transition. At finite chemical potential, G 2 
is determined by a flow equation resulting from the gen- 
eral requirement of integrability. By the flow of the ef- 
fective coupling, the basic features of the EoS at lx = 0, 
namely the nonperturbative behavior near confinement 
and the asymptotics, are mapped into the T-fx plane as 
exemplified by the Nf — 4 flavor system studied on the 
lattice. An important consequence of the quasiparticle 
approach is the relation of the critical values of temper- 
ature and chemical potential. For deconfinement matter 
with physical quark masses, this fact leads to the impli- 
cation that compact stars composed purely of /3-stable 
deconfined matter may be less massive and, hence, more 
distinct in the bulk properties to neutron stars than es- 
timated by other approaches. 

Acknowledgments: We are grateful to E. Crosse and 
F. Thielemann for initiating the present work. The stim- 
ulating interest of O. P. Pavlenko is acknowledged. The 
work is supported by BMBF 06DR829/1. 



[1] B. A. Freedman, L.D. McLerran, Phys. Rev. D16, 1130, 
1147, 1169 (1977) 

C. Zhai, B. Kastening, Phys. Rev. D52, 7232 (1995) 
G. Boyd, J. Engels, F. Karsch, E. Laermann, C. Leg- 
eland, M. Liitgemeier, B. Petersson, Nucl. Phys. B469, 
419 (1996) 

C. Bernard et al., Phys. Rev. D55, 6861 (1997) 
J. Engels, R. Joswig, F. Karsch, E. Laermann, M. 
Liitgemeier, B. Petersson, Phys. Lett. B396, 210 (1997) 
M. Le Bellac, Thermal Field Theory, Cambridge Univer- 
sity Press, Cambridge (1996) 

M. I. Gorenstein, S. N. Yang, Phys. Rev. D52, 5206 (1995) 
T. D. Lee, C. N. Yang, Phys. Rev 117. 22 (1960) 
A. Peshier, TFT'98 Proceedings, |hep-ph/9809379 
A. Peshier, B. Kampfer, G. Soil, O.P. Pavlenko, Euro- 
phys. Lett. 43, 381 (1998) 
E. Wang, U. Heinz, Phys. Rev. D53, 899 (1996) 

A. Peshier, B. Kampfer, O.P. Pavlenko, G. Soff, Phys. 
Rev. D54, 2399 (1996) 

B. Kampfer, O.P. Pavlenko, A. Peshier, M. Hentschel, 
G. Soff, J. Phys. G23, 2001c (1997) 
P. Levai, U. Heinz, Phys. Rev. C57, 1879 (1998) 
data files available from http://www.fz-rossendorf.de/ 
FWK/MITARB/Peshier 

[16] N. K. Glendenning, Compact stars, Springer Verlag, New 
York (1997) 



[4] 
[5] 

[6 

[7] 

[8] 

[9 

[10 

[11 
[12 

[13. 

[14 
[15 




FIG. 1. The characteristics of the coupling flow equation 
(H) for the QCD plasma with Nf — 4 light flavors for which 
G 2 (T, fi = 0) is obtained from lattice data |5J . At lead- 
ing order the characteristics are curves of constant coupling 
strength. The pressure is negative in the region below the 
dash-dotted line, thus excluding the region of intersecting 
characteristics. 
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FIG. 2. The effective coupling strength a = G 2 /(4-7r) as a 
function of fi and T for the Nf = 4 plasma in the chiral limit. 
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FIG. 3. Total pressure and energy density (lower and upper 
set of curves, respectively) of the charge-neutral, quark-gluon 
plasma in /^-equilibrium, scaled by the values of the free limit. 
The panels on the left (right) show the EoS at fj, = (T = 0), 
for values of the model parameter 3 < A < 9 (lower and upper 
line, respectively, of hatched area), and B = 120 MeV (top) 
and B — 180 MeV (bottom). Non-unique values of the 
energy only occur in the unphysical region, where p < 0. 



